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ABSTRACT 

We present optical spectroscopy of the microquasar SS 433 covering a significant fraction of a 
precessional cycle of its jet axis. The components of the prominent stationary Ha and H/3 lines 
are mainly identified as arising from three emitting regions: (i) a super-Eddington accretion 
disc wind, in the form of a broad component accounting for most of the mass loss from the 
system, (ii) a circumbinary disc of material that we presume is being excreted through the 
binary's L2 point, and (iii) the accretion disc itself as two remarkably persistent components. 
The accretion disc components move with a Keplerian velocity of > 600 kms -1 in the outer 
region of the disc. A direct result of this decomposition is the determination of the accretion 
disc size, whose outer radius attains ~8 Kq in the case of Keplerian orbits around a black 
hole mass of 10 M . We determine an upper limit for the accretion disc inner to outer radius 
ratio in SS 433, i?i n /-Rout ~ 0.2, independent of the mass of the compact object. The Balmer 
decrements, Ha/H/3, are extracted from the appropriate stationary emission lines for each 
component of the system. The physical parameters of the gaseous components are derived. 
The circumbinary ring decrement seems to be quite constant throughout precessional phase, 
implying a constant electron density of log 7V G (cm~ 3 ) ~ 11.5 for the circumbinary disc. The 
accretion disc wind shows a larger change in its decrements exhibiting a clear dependence on 
precessional phase, implying a sinusoid variation in its electron density log iV e (cm ) along 
our line-of-sight between 10 and 13. This dependence of density on direction suggests that 
the accretion disc wind is polloidal in nature. 
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1 INTRODUCTION 

Microquasars are X-ray binaries which undergo a wide range 
of physical processes including accretion onto a compact object 
(black hole or neutron star) and the launch of relativistic jets. 
SS433, one of the most studied microquasars, became famous as 
the first known source of relativistic jets in the Galaxy, and it is 
the only system, X-ray binary or active galactic nucleus (AGN), 
for which atomic emission lines have so far been associ ated with 
the jets, hence implying a baryonic content (i.e., e~p + , iMilgroml 
ll979tlCrampton & Hutchingslll 98 ll ; lender et all200d) . 

The optical spectrum of SS433 is characterised by the pres- 
ence of numerous broad emission lines with complex profiles, on 
top of a bright continuum. As at near-infrared and X-ray wave- 
lengths, these emission lines can be divided into two groups: lines 
that are referred to as stationary, albeit highly variable in strength 
and profile, and lines that are moving. The latter are thought to 
originate in the two oppositely-directed relativistic jets moving 
with a speed v ~ 0.26 c. The system SS433 shows four main 
periodicities: the binary's orbital motion, with a period of about 
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13.08 days dCrampton & Hutchings|[l98lh : the jet axis precession 
and n utation, with periods of about 162 and 6 days iKatz et all 
Jl982l) . respectively, a nd a recently discovered 550 day preces- 
sion of the radio ruff tooolin & Blund ell 2009). A configuration 
whe re two jets emerge in opposite directions from the central ob- 
ject l lFabian & Rees|[l979l) . in which the jet axis undergoes a pre- 
cession cycle every 1 62 days was pro posed to describe the motion 
of the emission lines jMilgrom|[T979h . and this is referred to as the 
kinematical model. 

The stationary optical and near-infrared spect rum of SS433 
is dominated by hydrogen and He I emission lines dMargon|[T984h 
and at least 15 per-cent of the flux i n the hydrogen lines is con - 
tributed by the accretion disc itself JPerez M. & Blundellll2009h . 
The emission-line spectra of accretion discs, and their evolution 
with orbital and precessional phases, comprise most of the infor- 
mation we can obtain about the temperature and density variation s 
as well as velocity gradients within the disc JSkidmore et al J2000h . 

The Balmer decrements of the stationary lines in SS 433 have 
not previously been studied due to the high interstellar extinction 
towards the object, which makes th e detection of the H/3 line rather 
difficult l lPanferov & Fabrikal 19971) . The decrements are highly de- 
pendent on physical parameters of the gas such as its tempera- 
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ture, optical depth and also the nature of the source of radiation 
dDrake&Ulrichlll98Cl) . 

The hard radiation field around compact objects and the ex- 
pected high electron den sities for SS433 (N e > 10 13 citT 3 , 
IPanferov & Fabrikd[l997h require that we use an adequate treat- 
ment for such environments in order to study the emission-line gas. 
At high densities, and in the presence of heavy ele ments, excitation 
by co ll isional processes becom e a relevant factor l lFerland & Reed 
1 1988b . iDrake & Ulrichl dl980t) performed theoretical calculations 
of the emission-line spectrum from a slab of hydrogen at moder- 
ate to high densities (10 8 < N c < 10 15 cm~ 3 ) over a wide range 
of physical para meters, including values c lose to those observed 
in SS433's gas l lPanferov & FabrikUl 19971) . We compare our esti- 
mates with Drake & Ulrich's findings in Section|4] 

The existence of dust mixed with the emitting gas would 
have e ffects on the Balmer decrem ent that are by no means neg- 
ligible dOsterbrock & Ferlanj2o"o3) . SS433's optical spectrum re- 
veals evidence of severe dust extinction, such as the remarkably 
red continuum and the presence of prominent i nterstellar abso rp- 
tion lines in the form of diffuse interstellar bands dMargonl 1984b . In 
our optical data described in sub sequent sections, and in agreement 
with previous ob servations (e.g. jMurdin et al .ll98d : lMargonl 19841 ; 
iGies et al]|2002l among others), we have detected the diffuse inter- 
stellar bands at 4430, 5778 and 5780 A and also the in terstellar lines 
Ca H A3968, Ca K A3934 as well as Na D A5890. iMurdin et al.l 
Jl980h reported an interstellar absorption (i.e., Galactic extinction) 
of Ay ~ 8 mag toward SS433, obtained fr om infrared meas ure- 
ments. This result has been corroborated bv lGiesetal.l (l2002h by 
fitting the spectral energy distribution, and the currently most ac- 
cepted value is Av = 7.8 mag. 

Fig. [T] shows the map of Galactic extinction in the direction 
of the SS433AV50 complex. We construct ed this extinction im - 
age by converting the colour excess map of ISchlegel et alj Jl998h 
into an estimate of the reddening, assuming a selective extinction, 
Rv = Av/E(B — V), equal to the a verage value for the Galaxy 
(i.e., R v = 3.1. ICardelli et alj|l989l) . where E{B - V) is the 
colour excess. It is clear from the image that the ~7.8 mag of 
extinction is consistent with the interstellar absorption gradient. 
The extinction estimate from this map towards SS433's position 
is Av = 7.813 mag. Although the value for E(B — V) is quite 
accurate, the choice of Rv is not. The intrinsic error introduced by 
Rv implies that the accuracy of the extinction estimate is not better 
than to a tenth of a magnitude. Therefore, we corrected our spec- 
tra (whose reduction is described in Section [2} for Ay — 7.8 mag 
of inte rstellar absorption using the reddening law of ICardelli et al .1 
dl989h . before calculating the Balmer decrements. 

In this paper we decompose the profiles of the stationary emis- 
sion lines, Ha and H/3, with a number of Gaussian components 
(Section |3)- Each model component is identified with its corre- 
sponding emitting origin. In Section|4]we analyse the physical con- 
ditions of each component present in SS 433 and its behaviour as a 
function of orbital and precessional phases, via the Balmer decre- 
ment. In Section [5] we present our calculations for the accretion 
disc radii and compare our findings with discs from other objects 
that share some common features with SS 433. 



2 OBSERVATIONS AND DATA REDUCTION 

We present optical spectroscopy data covering most of a preces- 
sional cycle of SS 433's jet-axis (and presumably therefore accre- 
tion disc). The data were acquired with the Supernova Integral Field 
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Right Ascension (degrees) 

Figure 1. Colour-scale image shows the Galactic extinction (Ay) to- 
wards the W50/SS 433 system from the IRAS/COBE all-sky survey 
lSchlegeletal.ll998l) . Colour scale corresponds to visual extinction in mag- 
nitudes and it is centred on SS 433's position (given by the red circle). Grey 
contours on the colour image are 3, 3.6, 5.4, 8.4, 12.6 and 18 mag. The 
radio continuum emissio n at 1465 MHz is shown as black contours (data 
from iDubner et all 1998). Radio contours are 12, 14 and 20 mjy beam -1 . 

Spectrograph (SNIFS) at the Univ ersity of Hawaii 2.2-m telescope 
by the Nearby Supernova Factory jAldering et al. 2002). The spec- 
trograph is composed of two modules, one for blue wavelengths 
covering the region between 320 nm to 560 nm with a resolution 
of ^1000 at 430 nm, and a module for red wavelengths (520 nm 
to 1000 nm) with a resolution of ~1300 at 760 nm. The data con- 
sist of 38 spectra (considering blue and red as one spectrum). The 
data set can be subdivided in two groups, one spanning from 2006 
April 27 to July 30 (precessional phases tp plc £ [0.7, 1.3]), and an- 
other one between 2006 October 1 and November 10 (precessional 
phases ip plc 6 [0.67, 0.95]). We use the convention in which or- 
bital phase ((ftor b) zero is when the d onor star is eclipsing out the 
compact object l lGoranskii et al.l 1998b . Prec essional phases are cal- 
culated based on the ephemeris reported in lFabrikal ( 120041) . where 
precession phase zero is when the moving jet lines are maximally 
separated hence the inclination of the jet axis with our line-of-sight 
attains a minimum, i.e., it corresponds to maximum exposure of the 
accretion disc to the observer. 

The data were reduced by the Nearby Supernova Factory 
group through their standard pipeline. Frequent observations of 
standard stars and the arc lamp were performed during each 
night to perform flux calibration and provide an accurate esti- 
mate of the wavelength axis. Wavelength calibration was carried 
out by fitting a high-order polynomial to the lamp spectra. Stan- 
dard star division was applied to each frame in order to flux cal- 
ibrate the data. Telluric features were not removed. All the sub- 
sequent analyses were carried out using the Perl Dat a Language 
(http : //pdl .perl . org, iGlazebrook et al.lll997l) . 

2.1 Fitting procedure 

Examples of the reduced spectra in the vicinity of the Ha and H/3 
lines are shown in Fig. [2] This figure also shows the best combi- 
nation of Gaussian components fitted to those spectra. The blue 
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Figure 2. Examples of H/3 (upper panel) and Ha (lower panel) so-called 
stationary emission lines. Both spectra were acquired on MJD 53009.23, 
which corresponds to ip pTe = 0.67 and </> or b = 0.69 (for a description of 
the convention used please see text in Section[2}. The different components 
are explained in the text. The narrow blue- (dot-dashed line) and red-shifted 
(dashed line) Gaussian components correspond to the accretion disc's flow. 
The top x-axis corresponds to rest velocity in units of kms -1 and both 
y-axes are in flux density units of W m~ 2 /im _1 . The broad wind compo- 
nent has a FWHM of ~ 1460 km s" 1 and ^lSOOkms^ 1 forH/3 and Ha, 
respectively. The error bars correspond to 3o\ which are about half the sym- 
bol size in the case of the Ha line. The shaded area represents that which 
we consider to be the wings of the accretion disc profile (see Section|5J. 

spectra were rebinned to the spectral resolution of the red spectra 
(130 kms -1 pixel -1 ). The underlying continuum was subtracted 
by fitting a low-order polynomial (linear in most cases and only in 
some cases parabolic). The best-fitting parameters were determined 
by x 2 minimisation in the pixel space. We explain the fitting model 
and the deconstruction of these profiles in the next section. 



3 DECOMPOSITION OF STATIONARY LINES 

The stationary emission lines in SS433 possess complex an d 
variable profi l es dBlundell et alj|2008l : IPerez M. & Blundellll2009t) . 
iFalomo et al.l dl987l) fitted the stationary Ha line profile of their 
low-resolution spectra using three to four Gaussian components. 
Their best fit consi sted of one broad co mponent plus three nar- 
rower components. iBlundell et alj d2008l) carried out decomposi- 
tion of the Ha line using higher-resolution data. They found one 
broad component (FWHM ~ 700 kms -1 ) whose width was ob- 
served to decrease with precessional phase (i.e., as the jets become 
more in the plane of the sky) identified as the accretion disc wind, 



and two narrower, red- and blue-shifted, components which were 
stationary in wavelength, being radiated from a glowing circumbi- 
nary ring. In the data presented in that paper, there was no evidence 
of the accretion disc lines being apparent, in contrast with their per- 
sistent appear ance in the data presented he re. 

Recently. IPerez M. & Blundelll J2009h were able to model the 
complex Brackett-7 (Br7) stationary line with usually 5 Gaussian 
components: a broad wind component (FWHM ~ 1500 kms -1 ), 
two Gaussian components accounting for the circumbinary ring 
and they also found the presence of two components being radi- 
ated by matter spiralling in the accretion disc. The presence of 
SS433's accretion disc in the near-infrared line seems to be per- 
sistent and their separation implies a circular velocity between 500 
and 700 kms -1 . We modelled our optical data following a similar 
approach. 

Following the analysis done by IPerez M. & Blundelll J2009h . 
the Ha and H/3 stationary lines were fitted with four Gaussian 
components in order to account for the complexity of the profiles. 
Although two narrow components have been used in the past to 
represent the split lines from the circumbinary ring dBlundell et al.l 
l2008l : |Perez M. & Blundellll2009h . because of the low resolution of 
the spectra we present in this paper (~ 1 30 km s~ 1 ) just one central 
narrow component was needed to represent the contribution from 
the circumbinary. Fig.[2]shows examples of fitted profiles. 

The presence of a P Cygni feature in the stationary lines 
at certain precessional phases has been noted by several authors 
dCrampton & Hutchingsl Il98ll ; iFilippenko et~al] 1 19881 ; iGies et al.l 
|2002|) . An absorption feature in the blue wing of the stationary 
line profile would indeed complicate the analysis and it would have 
to be taken into accoun t using models of outflowing winds (e.g., 
ICastor & Lamersll 19791) . We detected the P Cygni absorption sig- 
nature at the epochs ^ prc £ [0.7, 0.85] and we have excluded these 
data from our analysis. 

The reason why we identify the accretion disc with two sep- 
arate components representing each peak of the disc's profile is 
based on the model of emission lin e formation in accretion discs 
presented by iHorne & Marshl dl986l) . in which the velocity of the 
outer regions of the d isc is represented by h alf the velocity separa- 
tion of the line peaks dHorne &~M arsh 1986). This is caused by the 
Doppler effect and the rotational motion of the matter within the 
disc. Therefore, the speed with which the radiating material spirals 
within the accretion disc corresponds to half the difference between 
the speed at which those lines are moving (given by the relative po- 
sition of the peaks), under the assumption that the fitted centroids 
correspond to the tangent speed. This reveals material that is orbit- 
ing in the potential well at speeds of about > 600 kms" 1 (see blue 
and red tracks in Fig. [3}. 

Fig.[3]does not reveal a clear signature of the expected motion 
for th e accretion disc centro i d and certainly nothing corresponding 
to the iFabrika & Bvchkoval d 19901) sinusoidal plot that assumes a 
circular orbit. It is not possible to identify what would be the correct 
phasing of the motion of the disc due to lack of knowledge about 
the circularity or eccentricity of the orbit. For certain eccentrici- 
ties, there are not necessarily two well-defined peaks, for example. 
A detailed analysis of this problem requires much more long-term 
monitoring and sophisticated modelling. 

Fig. [3] (left and centre panels) clearly shows that the accre- 
tion disc lines appear slower at orbital phase zero. This is when 
the donor star obscures the inner (faster) region of the accretion 
disc. This is seen most clearly in Ha. Additionally, there appears 
to be an episode around MJD 53930, which corresponds to orbital 
phase 0.5, where the accretion disc lines seem to be rotating faster 
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Figure 3. Tracks of the fitted centroids for each component of the stationary H/3 (left) and Ha (middle) lines as a function of modified Julian day (j/-axis, 
increases vertically). The tick mark heights are proportional to the FWHM of each component. Red and blue ticks correspond to the accretion disc lines, the 
larger black ticks represent the disc wind while the light blue ones represent the circumbinary ring. The top x-axis corresponds to rest velocity in units of 
kms for an assumed systemic velocity of zero. The dot-dashed sinusoidal depicts the orbital motion of the binary. The dotted horizontal lines correspond 
to epochs when the orbital phase is zero. The plot on the right hand side corresponds to the evolution of the FWHM of the wind component of the stationary 
H/3 (gold triangles) and Ho (blue circles) lines. 



than usual and may simply correspond to a time when we have an 
unimpeded view into the inner, hence faster, part of the accretion 
disc. The right panel in Fig. [3] shows the evolution of the FWHM 
of the wind components of the stationary Ha and H/3 lines. This 
panel also shows that the wind lines appear to get broader at longer 
wavelengths, which is consistent with previous near-infrared ob- 
servations of the Br7 line, where the wind component has b een 
reported to be broader than in Ha jPerez M. & Blundellll2009l) . It 
is remarkable that near the end of the data set we present in this 
paper, the width of the accretion disc wind attains a velocity of al- 
most 2500 km s~\ which translates into a factor 1.7 increase in the 
10 -4 M ra yr 1 mass-loss rate estim ated from recombination line 
fluxes bv lPerez M. & Blundelll J2009I) . 



4 RESULTS 
4.1 Basic picture 



Modelling o f the stationary Balm er emission lines, Ha (this 
work and iBlundell et aU 120081) . H/3 (this work) and Br7 
jPerez M. & Blun dell 2009), yields a clear and consistent interpre- 
tation of the radiating components present in the SS433 system. 
The main components of the system, as inferred from the decom- 



position of the stationary line complexes, are: a rotating accretion 
disc, a fast outflowing wind that accounts for most of the mass loss 
in the system and a circumbinary ring probably being fed from the 
L2 point. The moving emission lines from the jets sometimes blend 
with the stationary ones at around ip ple ~ 0.3 or 0.4, contributing 
significantly to the intensity variation and profile shape. Moreover, 
care must be taken since the intensities and profiles of the s tationary 
lines vary strongly during flares (Blundell et al. in prep., iMargorj 
1 19841) . 



4.2 Balmer decrement of stationary lines 

The hydrogen Balmer lines show very variable Balmer decrements 
in SS433. Fig. [4] shows that the computed Balmer decrements for 
the stationary emission-lines scatter from 1 < Ha/H/3 < 8. In the 
int erpretation of th e se res ults we use the calculations carried out 
by brake & UlriclJ Jl980l) who modelled the emission-line spec- 
trum from a slab of hydrogen at high densities, by taking into ac- 
count important mechanisms such as collisional excitation and de- 
excitation, as well as self-absorption processes. It is important to 
chose reasonable temperature and optical depth (tl q ) parameter 
values in order to apply their calculations to SS433. We mainly 
utilise the decrements computed for a model with T c > 10 4 K. For a 
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Figure 4. The Balmer decrement, Ha/H/3, as a function of precessional 
phase (top rr-axis) or modified Julian day (bottom x-axis), for the accretion 
disc wind (top), circumbinary ring (middle panel) and accretion disc lines 
(bottom). The dot-dashed line corresponds to the best fitted sinusoidal to 
the disc wind decrements, with a fixed period of 162 days. In the bottom 
plot the red- and blue-shifted accretion disc lines are represented by solid 
triangles and circles, respectively. The vertical dotted lines represent times 
at which the orbital phase of the binary is zero. 



theoretically-modelled accretion flow, Kra mer's approximation fo r 
the optical depth yields tl^ = 4 x 10 torake&Ulrichlll980h . 
Since decrements for this exact value are no t reported in their paper , 
we use tl q ~ 10 5 instead. According to lDrake^fc Ulrichl J l980h 
these ratios (1 < Ha/H/3 < 8) then imply up to four orders of mag- 
nitude variation in density throughout the precession period. This 
large inferred variation is not a problem for plausible physical mod- 
els. We will analyse and discuss the Balmer decrements calculated 
for each component of SS 433 in the following sub-sections. 



4.2.1 Circumbinary disc decrements 

The Balmer line ratios of the componen ts correspondin g to the cir- 
cumbinary excretion disc observed by iBlundell et al. I d2008h arc 
shown in the middle panel of Fig. [4] They are roughly constant 
throughout precessional phase with an average value of Ha/H/3 = 
2.45 ± 0.69 (1 a scatter). This implies a moderate electron density 
of logTVc ~ 11.5, approximately constant throughout the preces- 
sional cycle and given the similarity of this Balmer decrement to the 
Case B canonical value of ~2.8, this is consistent with no attenua- 
tion from the accretion disc wind, i n accordance with the circ umbi- 
nary nature of this emitting region dPoolin & Blundellll2009l) . 



4.2.2 Accretion disc wind decrements 

Fig. [4] shows that in the case of the accretion disc wind, the 
Balmer decrements have a very clear precessional phase depen- 
dence, showing a clear tendency towards lower electron density 



(steeper decrements) as the jet axis becomes more in the plane of 
the sky. This dependence of density on direction is evidence for 
the polloidal nature of the wind. The decrement values for the disc 
wind vary between quite flat Ha/H/3 = 1.43 ± 0.28 between pre- 
cessional phases 0.95-0.05 and steep Ha/H/3 = 5.30 ± 1.75 for 
phases 0.65-0.84. These values correspond to variations in electron 
density of log N e ~ 13 and log iV ~ 10, respectively. 

4.2.3 Accretion disc decrements 

The accretion disc lines are the ones that show the most dra- 
matic and scattered fluctuations in decrement. Both lines, Ha and 
H/3, show a similar level of scatter. The average decrement for 
the red and blue components are Ha/H/3 = 2.78 ± 1.56 and 
Ha/H/3 = 3.34 ± 1.80, respectively. Inverted Balmer decrements 
can be seen at some epochs, implying very high densities and pos- 
sibly an enhanced inflow of matter from the companion onto the 
accretion disc. The decrement of the accretion disc itself appears 
to be dominated by fluctuations in the instantaneous wind density 
rather than being dominated by any precessional phase dependence. 



5 SIZE OF THE ACCRETION DISC 

A direct consequence of lHorne & Marshl 's l 19861 model of the emis- 
sion line profile formed in an accretion disc is that the emission 
close to the inner radius of the disc contributes to the outer wings of 
the double-peaked profile, while the regions neare r the outer radius 
provi de the emission for the peaks (see fig. 1 in iHorne & Marsh] 
1 1986b . By considering this along with the assumption that the ve- 
locity of the matter spiralling in the accretion disc i s Keplerian, we 
can estimate the size of the disc dMason et al . 2000). 

For a test particle in Keplerian motion at a distance r from a 
central object of mass Mbh, the Keplerian speed projected onto the 
plane of the sky is 



/ GMbh 
Vkep = \ I smt, 



(1) 



where i is the angle between the normal to the disc and our line- 
of-sight. The disc is edge-on when this inclination angle is 90° and 
pole-on when i = 0°. In our decomposition of the Balmer lines' 
profiles, the Keplerian velocity of the gas at the outer radius R out 
of the disc is given by half the separation between the peaks of the 
two Gaussian components representing the accretion disc, while 
the Keplerian velocity of the gas closer to the disc inner radius R ln 
corresponds to half the separation between the red wing of the red- 
shifted component and the blue wing of the blue-shifted component 
(see hatched area in Fig.|2j- 

The most suitable epoch to measure the inner and outer radii 
of the accretion disc is when the disc is in between the companion 
star and the observer (orbital phase 0.5) and when the disc is not 
totally edge-on (i.e., sini ~ 1) at precessional phase ip P i CC close 
to 0.33 or 0.66. The closest spectrum in our dataset to those epochs 
was taken on MJD 54009 which corresponds to ^ pre c = 0.67 and 
0orb = 0.69. This spectrum is depicted in Fig. [2] Both the half 
wing separation and the half peak separation for the accretion disc 
lines were measured from this spectrum. 

Half the separation between the peaks of the fitted Gaussian 
components of the disc yields a Keplerian velocity for the outer 
part of the disc of ~620 krns" 1 , measured from both the Ha and 
H/3 emission. The Keplerian velocity of the inner part of the disc 
was determined by measuring the positions in velocity of the line 
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Figure 5. Variation of the outer radius of the accretion disc as a function of 
the Keplerian velocity of the outer edge of the disc (half the separation of 
the line peaks, in km s — 1 ), for different values of the mass of the compact 
object. The solid area represents the range of parameters that we think are 
relevant to SS433. The thick orange line corresponds to iWbh = 10 Mq. 

wings (centre of the hatched areas in Fig. [2]l by inspection. Both 
Ha and H/3 present wing velocities of ~ 1150 kms - . 

The central compact object in SS433 is m ost likely to be a 
stellar mass black-hole, rather than a neutron star |Gies et al J2001 
lFabrikall2004 lLopez et alfeOOrj ; iBlundell et alj2008h . It is evident 
from Equation Q] that large black-hole masses imply larger outer 
radii for the accretion disc. The total mass of the accretion d i sc plus 
the compact object has been calculated by IBlundell et al. I d2008h 
and it attains 16 Mq. However, we point out that there has been a 
large dispersion re ported for the component masses of SS 433 (e.g., 
lLopez et al.ll2006h . We will assume that the mass of the black hole 
itself is of the order of 10 M©. 

The Keplerian velocity for the inner region of the disc of 
1200 kms -1 implies that the accretion disc's inner radius is 
R in ~ 1.4 M10 R©, where M w = M b h/10M Q . The estimated 
value for R[ n corresponds to an upper limit to the inner radius of 
the accretion disc at the time of observation. This is a direct conse- 
quence of the high obscuration presented in the system due to the 
accretion disc wind, which hinders the emission propagating from 
the inner parts of the disc. Additionally, the determination of the 
full-width at zero intensity (FWZI), or distance between the wings, 
could be imprecise sometimes and it could add certain degree of 
uncertainty to the R[ n value. Also, it is important to note that the 
inner radius of an accretion disc is known to recede as a conse- 
quence of changes in the accretion rate onto the compact object, 
for a range of masses of the compact object. 

On the other hand, the distance between the peaks of the ac- 
cretion disc components depicted in Fig.|2]yields a velocity for the 
disc outer region of about 620 kms -1 , which if Keplerian yields 
Rout ~ 5 M10 R©. Fig. [5] shows how Rout varies as a function 
of the Keplerian velocity of the outer edge. This is fully consis- 
tent with a previous estimate of the overall size of the system, 
whose semi-major axis of th e orbit has been reported to be ~80 R© 
JPerez M. & Blundellll2009h . 

The Keplerian velocity of the outer region of the accretion 

disc has been reported to be lower in the near-infrared Br7 line 
1 ; 1 11 -1 

dPerez M. & Blundell 2009) than in the optical lines analysed in 
this paper. This is consistent with observations of other accretion 
discs such as the one in WZ S ge. By analysing Doppler tomography 
maps, ISkidmore etafl d2000h found that the Br7 emission is being 



emitted from a region further away from the central object, imply- 
ing that this emission is more representative of the outer regions of 
the disc. This is probably the case in SS433 since the peak sepa- 
ration in Bt'7 is smaller than in the Balmer li nes. The reported half 
separa tion of the peaks of ~500 kms -1 bv lPerez M. & Blundelll 
d2009h implies an outer radius i? ou t = 8 M10 R© . 

The black-hole mass does not affect the inner to outer radius 
ratio, since Rin/Rout = (V ou t/Vin) 2 - In the case of SS433's disc, 
Rin/ Rout « 0.2 (for an outer radius i?. out ~ 8 M10 R©). 

The size of SS433's accreti on disc is qu i te larg e compared 
with other close binary systems. iMason et al. I ( 120001) determined 
the size of WZ Sge to be at most 0.43 R© for a white dwarf with the 
Chandrasekhar mass. This is about 20 times smaller than SS433's 
disc. The Ri n :R ou t ratio is approximately the same in both these 
systems. 



6 CONCLUSIONS 

We have presented optical spectroscopy of the microquasar SS 433 
covering a significant fraction of the precessional cycle of its 
accretion disc. The components of the prominent stationary Ha 
and H/3 lines have been identified as arising from three emit- 
ting regions: an accretion disc wind which is super-Eddington 
dPerez M. & Blun dell 2009), in the form of a broad component ac- 
counting for most of the mass loss in the system, a circumbinary 
disc of material probably excreted through the binary's L2 point, 
and the accretion disc itself as two persistent components, having 
an outer region Keplerian velocity of > 600 kms - . 

A direc t result of this decompositio n using our UKIRT data 
published in lPerez M. & Blundelll d2009h was the determination of 
the accretion disc size, whose outer radius attains 8 R©, for an 
assumed black hole mass of 10 M©. With the data presented in this 
paper we determined the accretion disc inner to outer radius ratio 
in SS433, i?i n /i?out to be ~0.2, independent of the mass of the 
compact object. 

The Balmer decrements, Ha /H/3, were extracted from the 
stationary emission lines for each component of the system. The 
decrement of the circumbinary ring seems to be quite constant 
throughout precessional phase, implying a fairly constant electron 
density of log iV c ~ 11.5 for the circumbinary disc. The accretion 
disc wind shows larger changes in its decrements as a function of 
precessional phase, implying variations in its density, N c , between 
10 10 and 10 13 cm . Thus, the physical parameters of the gaseous 
components imply rather dense environments emitting the Balmer 
lines. 



ACKNOWLEDGMENTS 

We thank the referee for a careful reading of the paper. We are very 
grateful to the Sciences and Technology Facilities Council (STFC) 
and Conicyt for the award of a STFC-Gemini studentship. We are 
specially grateful to Greg Aldering, the Nearby Supernova Fac- 
tory collaboration and the University of Hawaii for their generosity 
in allowing us to obtain some spectra on SS433 during intersti- 
tial time within their main supernovae follow up campaign. K. B. 
thanks the Royal Society for a University Research Fellowship. 

REFERENCES 

Aldering, G., et al. 2002, SPIE, 4836, 61 

© 2009 RAS, MNRAS 000,[T]r7] 



Blundell, K. M., Bowler, M. G., & Schmidtobreick, L. 2008, ApJ, 
678, L47 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 
245 

Castor, J. I., & Lamers, H. J. G. L. M. 1979, ApJS, 39, 481 

Crampton, D., & Hutchings, J. B. 1981, ApJ, 251, 604 

Doolin, S., & Blundell, K. M. 2009, ApJL, 698, L23 

Drake, S. A., & Ulrich, R. K. 1980, ApJS, 42, 351 

Dubner, G. M., Holdaway, M., Goss, W. M, & Mirabel, I. F. 1998, 

AJ, 116, 1842 
Fabian, A. C, & Rees, M. J. 1979, MNRAS, 187, 13P 
Fabrika, S. N., & Bychkova, L. V. 1990, A&A, 240, L5 
Fabrika, S. 2004, Astrophysics and Space Physics Reviews, 12, 1 
Falomo, R., Boksenberg, A., Tanzi, E. G., Tarenghi, M., & Treves, 

A. 1987, MNRAS, 224, 323 
Fender, R., Rayner, D., Norris, R., Sault, R. J., & Pooley, G. 2000, 

ApJL, 530, L29 
Ferland, G. J., & Rees, M. J. 1988, ApJ, 332, 141 
Filippenko, A. V., Romani, R. W., Sargent, W. L. W., & Blandford, 

R. D. 1988, AJ, 96, 242 
Gies, D. R., Huang, W., & McSwain, M. V. 2002, ApJ, 578, L67 
The Perl Data Language, Glazebrook K., Economou F., The Perl 

Journal, 1997, Issue 5, 5 - also reprinted in Dr Dobbs Journal, 

Fall 1997 issue 

Goranskii, V. P., Esipov, V. E, & Cherepashchuk, A. M. 1998, 

Astronomy Reports, 42, 209 
Home, K., & Marsh, T. R. 1986, MNRAS, 218, 761 
Katz, J. I., Anderson, S. E, Grandi, S. A., & Margon, B. 1982, 

ApJ, 260, 780 

Lopez, L. A., Marshall, H. L., Canizares, C. R., Schulz, N. S., & 

Kane, J. F. 2006, ApJ, 650, 338 
Margon, B. 1984, ARA&A, 22, 507 

Mason, E., Skidmore, W., Howell, S. B., Ciardi, D. R., Littlefair, 

S., & Dhillon, V. S. 2000, MNRAS, 318, 440 
Milgrom, M. 1979, A&A, 76, L3 

Murdin, P., Clark, D. H, & Martin, P. G. 1980, MNRAS, 193, 135 
Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of gaseous 
nebulae and active galactic nuclei, 2nd. ed. by D.E. Osterbrock 
and G.J. Ferland. Sausalito, CA: University Science Books, 
2006, 

Panferov, A. A., & Fabrika, S. N. 1997, Astronomy Reports, 41, 
506 

Perez M., S., & Blundell, K. M. 2009, MNRAS, 397, 849 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Skidmore, W., Mason, E., Howell, S. B., Ciardi, D. R., Littlefair, 
S„ & Dhillon, V. S. 2000, MNRAS, 318, 429 

This paper has been typeset from a TpX/ LTpX file prepared by the 
author. 



Balmer decrements in SS 433 7 



© 2009 RAS, MNRAS 000,[TJ{7] 



